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Abstract 

One-loop gap equations have recently been used by several authors to estimate the non- 
perturbative mass gap in a 3-dimensional gauge theory. I extend the method to two 
loops and demonstrate, that the resulting gap equation has a real and positive solution 
m — 0.34g 2 , which is in good agreement with the one-loop results and lattice data. 
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In the high-temperature phase of the Standard Model one naively expects a vanishing Higgs 
vacuum expectation value and vector boson mass. This leads to the well-known breakdown of 
the perturbative expansion due to severe infrared divergences in the magnetic sector of the 
theory The infrared problem may be cured by a non- vanishing magnetic mass, which 
acts as a cutoff regularizing these divergences. Its inverse, the magnetic screening length, 
determines the size of non-perturbative effects in the symmetric phase and is closely related to 
the confinement scale of the effective 3-dimensional theory which describes the high-temperature 
limit of the 4-dimensional finite temperature field theory [[|. In an apparently massless 3- 
dimensional Yang-Mills theory the gauge coupling g 2 carries the dimension of mass, thereby 
providing a natural mass scale. A popular framework for calculating the dynamically generated 
mass in a 3-dimensional SU(2) gauge theory has recently been via gap equations. 

Up to now, several attempts have been made to evaluate the size of this mass gap, all of 
them at the one-loop level [[^,|],^]||. In order to check, whether the whole approach is consistent, 
it is crucial to extend the method to two loops. To my knowledge this letter contains the first 
treatment of gap equations in next-to-leading order. The corresponding calculation is rather 
involved at the technical level. 

After briefly presenting the idea behind the gap equation approach and giving an overview 
of the one-loop results, I will give reasons for the necessity of a two-loop calculation, the results 
of which will be extensively discussed. 

General Strategy 

All the models considered so far start form a massless Yang-Mills action Sg in 3-dimensional 
Euclidean space and then add and subtract some gauge- invariant mass term S m , where the 
subtracted term enters the perturbation theory at one loop higher than the added term. This 
can be formalized by introducing a loop-counting parameter / in the following manner: one 
rescales all the fields by y/l and calculates with the modified action 

Seff = y (S G (VlW) + S m {VlW)) - S m {VlW) (1) 

in a formal /-expansion || . Perturbative calculations are no longer done to a fixed order of the 
gauge coupling g, but as a power series in /, resulting in a rearranged or resummed perturbation 
series. Out of the many choices, which are possible for S m , I will concentrate on the non- linear 
cj-model for reasons specified below. 

The gap equation is a self-consistent condition for the vector boson mass. The goal is to find 
the particular size of the tree- level mass-term m = Cg 2 leading to a convergent perturbation 
series. In other words, the pole of the transverse part of the (Euclidean) vector boson propagator 
should remain at p 2 = —m 2 to any loop-order, i.e. 
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with some residue Z. With eq. (@) one obtains the desired gap equation for the self-energy in 
resummed perturbation theory 

IW = -m 2 ) (l + ^(p 2 = -m 2 )) = . (3) 

In n-th order of resummed perturbation theory one calculates eq. (^) up to l n and solves the 
gap equation for m. 

At one-loop eq. (Q) reduces to 

U l-loo P(j) 2 = _ m 2) _ g (4) 

In theories with a BRS-symmetry the position of the pole of the propagator and therefore 
eq. (|3]) is gauge- independent on mass-shell [[?[]. The self-energy itself is not gauge-invariant on 
mass-shell except at one-loop level. 

One-Loop Calculations 

In recent years, several authors have proposed models to extract a gap mass at one-loop 
level. When studying the electroweak phase transition Buchmuller and Philipsen could obtain 
a non- vanishing vector boson mass in the symmetric phase of a linear 3-dimensional SU(2)- 
Higgs model 0. A mass resummation was supplemented by a vertex resummation in order to 
get a BRS-invariant resummed tree-level action resulting in a gauge-independent gap equation. 
Deeply in the symmetric phase the value for the gap mass is approximately the same as the 
one obtained in a non-linear a-model. This suggests to investigate first the two-loop effects in 
this simpler model. 

In order to minimize the amount of diagrams in a two-loop calculation I follow a suggestion 
by Jackiw and Pi 0. The functional integral for the partition function in the non-linear a- 
model, where S m = S a , reads 

Z = J DWDnAexp - -(S G + S a + S GF - lS a ) , (5) 

with Sqf being some gauge-fixing term, which depends only on W M , and A the corresponding 
Fadeev-Popov determinant. They integrated out the Goldstone and ghost fields exactly (in 
an arbitrary gauge) and arrived at a massive Yang-Mills theory without any additional gauge 
fixing terms, 



Z oc J DWexp -- ^S G (V~IW) + m 2 Tr J d^VlW^VlW^ - lm 2 Tr J (f xVlW^VlW^ 



(6) 



One can also view this massive Yang-Mills theory as a non-linear cr-model (with i?£-gauge fixing) 
in unitary gauge. A calculation of the self-energy to one-loop in both theories indeed yields the 
same result on mass-shell and therefore the same gap mass. The off-shell self-energies coincide 
only in the limit £ —>■ oo. In unitary gauge the result for the off-shell transverse self-energy is 



l07T 



p 6 2p A 10p 2 \ 2m p p 4 Ap 2 
+ 8 arctan 1 



4m 6 m 4 m 2 I p 2m m 4 m 2 



(7) 



The gap equation for the massive Yang-Mills theory to one-loop is then a linear equation 
for m and reads 

1 63 

^m 2 SM - -—(—\n3 - 3) lg 2 m SM = , (8) 

1071" 4 

or 

m SM ^ 0.28/ . (9) 

Note that in unitary gauge the longitudinal part of the self-energy vanishes for all external 
momenta, n L (p 2 ) = 0. 

For the 3-dimensional SU (2) gauge theory also other gap equations have been considered, 
which are based on the Chern-Simons eikonal and on the non-local action 

= m 2 Tr J d 3 xF^F, , (10) 

where F p = \e pvp F up Interestingly, the one-loop gap equation of Alexanian and Nair 

yields a magnetic mass closely related to msM, 

4 2 
m A N = ^ m SM - 0.38# . (11) 

Jackiw and Pi obtain a complex magnetic mass with the mass term of eq. (|T0|) , which, however, 
can be modified such that the generated mass gap becomes real. Another attempt was recently 
made by Cornwall 0. His pinch-technique gap equation led to a mass gap of 

m c ^ 0.25g 2 . (12) 

It is also very encouraging, that these analytically calculated gap masses are consistent with 
the propagator mass obtained in a numerical lattice simulation in Landau gauge ||, 

m L SM = 0.35(1)/ . (13) 



Two-Loop Gap Equation 

Why is it crucial to perform a two-loop calculation? 

• The loop expansion does not correspond to an expansion in a small parameter. Neverthe- 
less, it might very well be, that the one-loop results provide reasonable approximations 
of the true mass gap. This can only be clarified by a two-loop calculation. 



• If the whole method is consistent, the numerical values for the mass gap in the different 
models should converge at higher loop-orders, since to all orders they describe the same 
Yang-Mills theory. In this case one expects the two-loop correction to be of order m^jy — 

• The two-loop gap equation is quadratic in m, whereas at one loop it is linear. The 
existence of a positive solution is a non-trivial check of the whole approach. 

As mentioned above, we first consider the Lagrangian of a massive YM (a resummed non-linear 
a-model in unitary gauge) 



C = -F* Ft, + -m 2 W°:W°: - -m 2 WyV a a , (14) 



\KvKv + -m 2 W a a W a a - 
with 

F^ = d,W^-d u W; + Vlge abc W b ^, G = SU(2), d = 3 - 2e . (15) 

Eq. (^) has to be expanded up to 0{l 2 ), which requires the evaluation of the diagrams 
depicted in fig. |I|. The contribution of diagrams 1, 2 and 3 to the transverse on-shell self-energy 
was already calculated at one loop, cf. eq. (|8|). The contributions of diagrams 4 and 5 are easily 
evaluated, 

tf-loop-CT^ = _ m 2 } = _L ^ ln3 _ g^j p g 2 m (16) 

For T^Ilr, which contributes to 0(1^), only one-loop diagrams are needed in the two-loop gap 
equation, which can directly be obtained from eq. (0), 

^IT^ 00 V = -m 2 ) = -L (33 - -ln3) I 9 - . (17) 
dp 2 T KF ' 32tt V 4 J m K J 

Far more work has to be done for the evaluation of the remaining 9 two-loop diagrams, 
which contribute to 0(l 2 g 4 ). As all propagators are massive and the external momentum 
does not vanish, the reduction of the scalar integrals to basic integrals with no momenta in 
the numerators turns out to be the most difficult step in the calculation. For propagator 
type integrals this task has been achieved only recently by Tarasov fLOf . Using his recurrence 
relations it is possible to reduce the self-energy integrals to a small set of linearly independent 
basic integrals. For the first time this method achieves a complete reduction and stays on an 
algebraic level as far as possible. Since the recurrence relations are in some cases quite involved, 
they have to be implemented into a FORM package In unitary gauge the situation is even 



more complex due to the high powers of momenta in the numerator. 

In 3 — 2e dimensions, the result of the reduction to basic integrals for the on-shell self-energy 
in unitary gauge reads 
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Except for — Q) - > which has to be evaluated numerically, there exist analytic expressions 
for the basic integrals in 3 — 2e dimensions JT2] ]. The result of ( |T8"D was also obtained using 
a FORM package written independently by O. Tarasov. Neglecting the resummation counter- 
terms, which spoil BRS-invariance, I also calculated the self-energy for the non-linear a-model 
in Feynman gauge, £ = 1. I obtained the same position for the pole of the propagator as in 
unitary gauge, which constitutes a very stringent test for the algorithm I used. 

Two further remarks have to be made concerning the two- loop calculation. First, the 
longitudinal part of the self-energy in unitary gauge vanishes for all momenta at two-loop level, 
which is another nice check of the calculation. Second, the non-linear sigma model is non- 
renormalizable. This is not a problem at one loop, since the 3-dimensional self energy is finite 
in dimensional regularization. At two loops, however, CD and Q are UV- divergent, 
which requires the addition of counter-terms (~ I 2 ) to the Lagrangian. The explicit calculation 
in Feynman gauge shows that a mass and wave function renormalization is sufficient to remove 
the infinities in the self-energy, 

J-Ul-'lf{p 2 ) = (L-Jl^-) - JL + finite . (19) 
ZV T ' c_1 V ' \12 60m 2 / 647r 2 e V ' 

Compared to the unitary gauge, where ghosts and Goldstone bosons are integrated out, the 
Feynman gauge involves the evaluation of many more diagrams, 33 generic two-loop graphs 
instead of 9. Note also, that the unitary gauge is not suitable for renormalization. Even 
in renormalizable theories the bad high-energy behavior of the propagator leads to terms 
(m 5 ) e' n > 1> m ^ ne se lf- ener gy> which cannot be dealt with by a mass or wave function 
renormalization. Renormalization in Feynman gauge introduces a renormalization scale //. Us- 
ing the MS'-scheme it turns out that there is almost no numerical dependence of the two-loop 
gap mass on the scale \i = fJ-j^g. 

Now we are ready to discuss the two-loop gap equation. In an arbitrary gauge of the 
resummed non-linear a-model it reads 



lm 2 - 0.28455 lg 2 m + f^) l 2 g 2 m + / 2 (£) l 2 g 2 m 

- 0.064346 ZV + 0.0037995 l 2 g A ln-^ = , (20) 

m 



with 
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The two- loop gap equation is not exactly gauge parameter independent. There remains a 
weak gauge dependence stemming from diagrams which involve the resummation counter-terms 
(mass counter-terms for the vector, ghost and Goldstone field). In the unitary gauge f\ and 
/2 reduce to eq. (16) and Note that the limit £ — > oo has to be taken before divergent 



integrals are evaluated The gap equations in unitary and Feynman gauge turn out to be 
identical. 



Two-Loop Results 



The major result is that the gap equation (|20|) has indeed a real and positive solution for 
p. The results are given in table 1 for different values of /x and £. The two-loop correction to 
the one-loop gap mass (|9|) is only 15 — 20%. 



JL 
m 


0.3 


1 


3 




0.343 


0.335 


0.327 


3 ' s 


0.345 


0.336 


0.328 




0.350 


0.342 


0.334 



Table 1: Solutions of the two-loop gap equation 

One may worry about the dependence of the gap mass on the renormalization scale \x and 
on the gauge parameter £. This is an artefact of (resummed) perturbation theory, which is 
expected to be cancelled at higher orders. Fortunately, the dependence of ^ on /i and £ is 
numerically unimportant. This suggests that the solution constitutes a reliable approximation 
to the exact gluon propagator mass in SU (2) gauge theory. 

A vector boson mass ~ 0.34g 2 is not in contradiction with confinement. It is of the same 
size as the confinement scale given by the string tension which was calculated in [Q. The 
connection of such a propagator mass to the heavier glueball masses ~ 0(l)g 2 [|T3| in a 3- 
dimensional SU(2) gauge theory remains to be clarified. 

The two-loop result survives all the crucial tests which have been mentioned above unex- 
pectedly well: the quadratic gap equation has a real and positive solution, which is not far 
away from the one-loop result. Moreover, it is now in better agreement with Alexanian and 
Nair's gap mass and matches perfectly the lattice result obtained by Karsch et. al. 

To judge the significance of a calculation in the non-renormalizable non-linear sigma model, 
it is crucial to perform the whole calculation in the linear Higgs model, which is super- 
renormalizable. The linear model involves the summation of hundreds of two-loop diagrams. I 
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have shown that the pole of the Higgs and of the vector boson propagator is gauge-invariant to 
two loops, which is a very powerful test of Tarasov's algorithm and my FORM package. More 
importantly, the two-loop gap equation in the linear Higgs model turned out to be numerically 
nearly independent of the Higgs mass. Therefore, the non-linear sigma model remains a very 
good approximation, as it already proved to be the case at one loop. Detailed results of this 



work will be published in a forthcoming paper 16 



I thank O. Tarasov for an independent calculation of eq. QTSI) as well as for helpful dis- 
cussions. I am also grateful to W. Buchmuller and O. Philipsen for valuable comments and 
suggestions. 
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